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Phase/gain imbalance estim ation or compensation 



5 

The present invention relates to a receiver for estimation or compensation 
of phase imbalance or gain imbalance and methods which can be applied in a receiver 
for estimation or compensation of phase imbalance or gain imbalance. A transmission 
technique used in the receiver is based on a QPSK (quadrature phase shift keying) 
10 modulation and a modulation scheme based on complex scrambling code. 

The Universal Mobile Telecommunications System or UMTS is one of 
the major new 'third generation 1 (3G) mobile communications systems being developed 
within the framework defined by the ITU (International Telecommunication Union) and 
known as IMT-2000. 

15 

Methods for estimation or compensation of phase imbalance or gain 
imbalance typically rely on reference signal generation, adaptive interference 
cancellation, blind source separation or a least square approach in the Fourier domain. 
Such methods are known, e.g. from "Digital I/Q imbalance compensation in a low-if 

20 receiver", Jack P. Glas (Bell Labs)-1998-IEEE, 'Improving QPSK Demodulator 
Performance for Quadrature Receiver with Information from Amplitude and Phase 
Imbalance Correction", Hung Ngyen (The Aerospace Corp.)-2000-IEEE, "Blind Source 
Separation Based I/Q imbalance compensation", Valkama, Renfors, Koivunen-2000- 
IEEE, "Advanced methods for I/Q imbalance compensation in Communication 

25 Receivers", M, Valkama, M.Renfros, V.Koivunen-DBEE transactions on signal 
processing- Vol. 49, no. 10, Oct 2001, 

WO 0150616 discloses a receiver for correcting the phase error in a 
received signal. However, this receiver does not provide an efficient estimation or 
compensation of phase imbalance or gain imbalance in the received signal. 



30 
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It is an object of the present invention to provide an efficient estimation 
or compensation of the phase imbalance or the gain imbalance. 

According to an exemplary embodiment of the present invention, this 
object is solved with a receiver for estimation or compensation of phase imbalance or 

5 gain imbalance utilizing a QPSK modulation scheme and a modulation scheme based 
on complex scrambling code, the demodulator comprising means for estimating the 
phase imbalance or gain imbalance before synchronisation. 

Such a receiver according to an exemplary embodiment of the present 
invention as set forth in claim 1, allows the estimation or compensation of phase 

1 0 imbalance or gain imbalance that can affect the overall performance of the receiver, at 
an early stage, i.e. before synchronisation, by e.g. exploiting the properties of the 
complex scrambling code. Thus, phase imbalances or amplitude imbalances will not be 
able to introduce losses in the further phases of the connection. Furthermore, if the 
bandwidth inside which time variations in the imbalances occur is lower than the 

15 transmission bandwidth, the demodulator is also able to track these time variations. 

The phase/gain imbalance estimation and compensation according to the 
present invention is e.g. tailored to the UMTS modulation scheme, i.e. may exploit the 
special properties of the UMTS modulation scheme for estimation and compensation of 
imperfections like phase-gain imbalances in the RF receiver or demodulator in an early 

20 stage of the receiver. 

According to another exemplary embodiment of the present invention as 
set forth in claim 2, 1 and Q components of the received signal signal affected by phase 
imbalance or gain imbalance after demodulation are determined. A ratio between a 
cross correlation of the I and Q components (<IQ>) and the mean value of the square of 

25 the I component (<P>) or a ratio between a cross correlation of the I and Q components 
and a square root of the product between the mean value of the square of the I 
component and the mean value of the square of the Q component ((<3 2 ><Q 2 >) 1/2 ) or a 
ratio between the mean value of the square of the Q component, (<Q 2 >) and the mean 
value of the square of the I (<P>) component are determined. 

30 According to yet another exemplary embodiment of the present invention 

as set forth in claim 3, a low pass for low pass filtering the signals are being 



implemented. By doing so, errors in the estimation of phase imbalance or gain 
imbalance can very efficiently be reduced, while computing mean values. 

According to yet another exemplary embodiment of the present invention 
as set forth in claim 4, the receiver comprises means for compensating the phase 
5 imbalance or gain imbalance before synchronisation. 

According to another exemplary embodiment as set forth in claim 5, the 
receiver is a WCDMA (UMTS) receiver (Wide Band Code Division Multiple Access 
receiver). Thus, as a basic principle exploited in the WCDMA (UMTS) receiver, the 
channel is frequency selective with respect to the transmission bandwidth allowing for a 
10 very accurate estimation of the phase/gain imbalance before synchronisation. 

According to another exemplary embodiment as set forth in claim 6, the 
estimation of the phase imbalance or gain imbalance is carried out iteratively providing 
for very accurate results. 

According to yet another exemplary embodiment as set forth in claim 7, a 
1 5 method for estimation or compensation of phase imbalance or gain imbalance in a 
receiver is provided, which may advantageously exploit the properties of the complex 
scrambling code and of the overall UMTS modulation scheme in the early stage of the 
receiver by estimating the phase imbalance or gain imbalance before synchronisation. 

According to an exemplary embodiment of the present invention as set 
20 forth in claim 8, a very simple method is provided, requiring only minimized 

computation power, and which may advantageously allow an iterative solution for 
estimating compensating both phase and gain. 

According to yet another exemplary embodiment of the present invention 
as set forth in claim 9, a feed-forward scheme or a feed-back scheme is established. 
25 According to another exemplary embodiment as set forth in claim 10, the 

estimation of the phase imbalance or gain imbalance is carried out iteratively providing 
for a very efficient compensation/estimation. 

According to the present invention, there is also provided a computer 
program for estimation or compensation of phase imbalance or gain imbalance as set 
30 forth in claim 11. 

It may be seen as a gist of the present invention that the phase imbalance 
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or gain imbalance is estimated or compensated before symbol synchronisation. Thus the 
phase imbalance and gain imbalance will not be able to introduce losses in the further 
phases of the connection; also the method is able to track temporal drifts in phase and in 
gain if the bandwidth inside which these phenomena occur is lower than the 
5 transmission bandwidth. The advantage of this method is, that it can be used both, for 
phase imbalance estimation and compensation and gain imbalance estimation and 
compensation in the case that there is also a gain mismatch. 

These and other aspects of the present invention will be apparent from 
and elucidated with reference to the embodiments described hereinafter. These 
1 0 embodiments will be described with reference to the following Figures: 



Fig. 1 shows a feed-forward scheme for phase imbalance estimation and 
compensation according to an aspect of the present invention. 
15 Fig. 2 shows a feed-back scheme for phase imbalance estimation and/or 

compensation according to an aspect of the present invention. 

Fig. 3 shows a feed-back scheme for gain compensation according to an 
aspect of the present invention. 

Fig. 4 shows schematically preferred application points for the method 
20 according to the present invention. 

Fig. 5 shows simulation results for different filter lengths and a velocity 

of3km/h. 

Fig. 6 shows simulation results for different filter lengths and a velocity 

ofl20km/h. 

25 Fig. 7 shows simulation results for different filter lengths and a velocity 

of 250 km/h. 

Fig. 8 shows the effect of the phase imbalance on the BER. 
Fig. 9 shows simulation results for the square root error for 5 degrees 
initial phase error vs. low pass integration length. 
30 Fig. 10 shows an estimation for a time varying phase (5 Hz). 
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Fig. 1 1 shows an estimated phase in accordance with an iterative 
implementation of the method according to the present invention. 

Fig. 12 shows a gain error in accordance with an iterative implementation 
of the method according to the present invention. 



5 



In the following description of the above Figures, the same reference 
numbers are used for the same or corresponding elements. 

Figure 1 shows a device with a feed-forward scheme for phase 
imbalance estimation and compensation according to an aspect of the present invention. 
As will be shown below, the demodulation method according to which the device of 
Fig. 1 is operating relies on the fact that the correlation matrix of the transmitted 
components, (i,q) is diagonal. As consequence, if I and Q are the in-phase and 
quadrature components of the received signal (with I = k • i and Q = k • (-isincp + 
q-coscp) as shown in equations (1)), the correlation between I and Q, even after 
propagation through the fading multipath scenario, comprises only a term due to phase 
imbalance. This is due to the property of the QPSK modulation scheme and the 
complex scrambling code. This is shown in more detail in the following with respect to 
a the UMTS system: 

The waveform transmitted in accordance with UMTS 
The UMTS transmission scheme in the downlink is based on a: 



transmitted to a symbol of a quaternary constellation. QPSK symbols can be tough as 
complex numbers, and the application of the complex scrambling code to the QPSK 
symbols represented in terms of multiplication between complex numbers. Ih other 
words, to scramble the symbol (at the instant n) s(n)=<sl(n),s2(n)), the symbol for the 
scrambling code element c(n)=(cl(n),c2(n)) has to be multiplied in accordance with the 
following rule: 



QPSK modulation 
Complex scrambling code 

The QPSK modulation scheme basically maps each two chips to be 



s • c s (sl(n) • cl(/i) - s2(n) • c2(n), sl(n) • c2(n) + s2(n) • cl(/i)) 



(1) 
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The two sequences composing the complex scrambling code, i.e. the 



5 



sl(n) and s2(n) are selected in a family of Gold codes and have the following properties 
in first approximation: 

sl{n) * sl(n + /w)=0,m*0 

sl(n)*s2(n + m) = 0 9 Vm (2) 
where * refers to the correlation operation. 

It is important to note that above equations (2) hold true in first 
approximation, and that the correlation interval can be extended to the full length of the 
code, but also to sub-intervals. 

Indeed, if the components i and q of the transmitted signal given from the 
equation (1) are considered, the following equations (3) can be derived: 



sequences, but also for intervals included if the included intervals are selected to have 
sufficient length. In the case of the UMTS, for the bandwidth involved, it is common 
understanding that the sufficient length can be extended to about 256 chips, i.e. the 
duration of a pilot symbol. 

The in-phase and the quadrature components of the received signal are 

uncorrelated 

In the following, it is shown that after propagation through a multi-path 
fading environment, if there is no phase imbalance, the in-phase and the quadrature 
components of the received signal are still uncorrelated. 

At this point, i.e. before decoding, the signal is a mixture of voice, date 
and control information. The signal at the receiver can be written: 



where the channel attenuation, delay and phase rotation is considered for each path. In 
what follows, I and Q are the in-phase and quadrature component of the received signal. 
Ii and Qi are the received components on the i-path, and i, (t - ) and q t {$ - r g ) 
represent the i-esima replica of the in-phase and quadrature components of the 



i(n) = sl(n) • cl(n) - s2(ri) • c2(n) 
q(n) = sl(n) • c2(n) + s2(n) • c\(n) 
i(n) * q(n) =< i(n) • q(n + m) >= 0 



(3) 



These relationships hold true for the maximal length of the Gold 



K0 = 2>,-e M *('-*,) = Z'i(') 
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transmitted signals i(t) and q(t), due to the path i. r f is the delay of the transmission 
path. 

The I and Q components of rj(t) can be expressed: 
Ii = &[ii(t-Ti>cos(q>i) + qi(t-Ti>sin(q>i)] 
5 Qi = gi[ii(t-Ti)(-sin(90) + qi(t-Ti>cos(cpi)] 

from which it follows for the elements of the auto correlation matrix : 
<L?> = <© 2 > • <i?> ^cosCcpO^ + <gi 2 > • <qi 2 > ^sinfaO^ - 
2<gi 2 > ■ <iiqi> •<sin(<pi)cos(<p i )> 

<Qi 2 > = <gi 2 ><ii 2 > -<sin(cpi) 2 > + <gi 2 > • <qi 2 > ^cosfaOS* - 
10 2<g 2 > • <iiqi> <sin((pi)cos(<pi)> 

<Ji Qi> = <g 2 > • -<-sin((pi) cos(9i)> + <g 2 > • <qi 2 > -<sin(90 cos(9i)> - 
<g 2 > • <iiqj> -<sin((pi) 2 > - <gi 2 > • <iiqi> -<cos(cp0 2 > 

It has to be noted that <ij 2 > = <qi 2 > = P/2 holds true for each I, with P 
being the transmitted power for each symbol. 
1 5 Also, <ii-qi> = 0 holds true for the orthogonal property of the in-phase 

and quadrature components of the transmitted signal . 

Also, the terms g, representing the channel attenuation for each path i, 
depend on the time. However, their dependence is slow with respect to the period onto 
which the mean operations are extended, and it can be assumed that they are 
20 uncorrected with the temporal variations in the transmitted signal. As consequence, the 
channel terms in the mean operations can be separated and written as follows: 
<& 2 >=gi 2 (t) 

The same separation is performed for sin(9i) and cos(9i) giving the phase 
information for each path. 
25 It follows that 

<I i 2 > = <Q i 2 > = gi 2 (t)*P/2 

i.e. the received power is determined by the transmitted power and the channel gain. 
Furthermore, it follows that 

<J r Qi> = 0 

30. i.e. the signal received on one path has in-phase an quadrature components that are 
orthogonal. 
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Assuming a normalized power (per bit) during transmission (i.e. P/2=l), 
a consideration of the correlation between in-phase and quadrature components on the 
same path yields the following matrix: 

^Qr'Ir) <Qr 2 > J U ^(0, 

where: 

The correlation between the terms due to different paths 
In the following, the correlations between the terms due to different paths 
are considered. 

The cross correlation between elements of different paths can be written 

as follows: 

<Ii-Ij> = <grgj> • <ii,ij> -<cos(<pi) cos(q>j)> + <gigj> • <qiqj> -<sin(9i) sin(q)j> - 
<gi,gj> • <iiqj><sin(<pj) cos(q>0> + <gigj> • <q i i j xcos(9 i ) sin(9j)> 
<Ii- Qj> = <gigj> • <ii4j> <-sin(cpi) cos(q)j)> + <gjgj> • <qjq,> <sin(9i) cos(9j> - 
<gi -gj> • <irqj><sin(9i) sin(9j)> + <gi-gj> • <qrijXcos(9i) • cos(9j)> 

where, in the temporal mean, the terms due to the channel, gi and 9^ and 
the statistical operations involving these terms have been separated. Similar expressions 
can be derived for <QrQj> and <Qj-Ii>. 
Since: 

<ij-ij> = <qrqj> = <qi-ij> = 0 

for i 56 j, all these terms are zero both in the short term and in the long term. 

Therefore it can be concluded that the correlation matrix between the 
components I and Q at the receiver is diagonal. Furthermore, elements on the principal 
diagonal depend on the time varying channel. However, it can be derived from equation 
(7) stated below, and from the fact that the channel is varying slowly with respect to the 
temporal extensions on which the temporal means for computing the correlations are 
taken, that this dependency does not influence the estimation of the gain/phase 
imbalance, since it is present at both, the numerator and at the denominator of the final 
expression for sin(9). 
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As shown in Fig, 1, in a feed-forward scheme for phase imbalance 
estimation and compensation according to an aspect of the present invention, an 
incoming complex signal r(t) is received at a receiver station not shown in Fig. 1 and is 
further processed such that the complex signal r(t) is provided to input connections 1 
5 and 2 of the receiver or demodulator 3. After that, the complex signal r(t) is 

demodulated. Between the input connections 1 and 2, the signal A cos(cot) represents 
the waveform used in the receiver to demodulate the incoming signal. Signal I on line 4 
and signal Q on line 5 are called the I and Q components of the demodulated complex 
signal r(t) in the following The I signal or component is the "in-phase" component of 

10 the signal r(t) and the Q signal or component is the quadrature component. 

The I and Q components affected by phase imbalance or gain imbalance 
and are provided to lines 6 and 7, which are connected to lines 4 and 5. The I and Q 
components are fed into means 8 for taking the cross correlation of I and Q, such as a 
multiplier. Additionally, the I component is fed into means 9 for taking the quadrature 

15 of I, such as a multiplier which multiplies I with I. After that, the resulting signals are 
taken through low pass filters 12, 13 by lines 10, 11, respectively. For estimation of the 
phase imbalance, a ratio between the cross correlation of I and Q and the mean value of 
the square of I is generated by means for generating a ratio 14, such as a divider, and 
consecutively passed through low pass filter 15. After that, the resulting signal is 

20 provided to a phase/gain compensator 17 via a line 16. The phase/gain compensator 17 
compensates the phase imbalance or the gain imbalance of the I and Q components. 
Thereupon, the estimated and compensated I and Q components i* and q c , are taken by 
lines 18 and 19, respectively, to be fed into means for synchronisation 20, such as a 
UMTS synchronizer. After synchronisation, the signal is ready for further processing, 

25 e.g., in reference to the standard architecture for the UMTS receiver, it can be fed into 
the rake receiver to exploit time diversity and then to the channel decoder to 
correct/reveal residual errors. 

Reference character 21 designates a multiplier or squaring device 
connected to line 5. An output of the multiplier 21 is connected to a low pass 22. An 

30 output of the low pass 22 is connected to another divider 23. The divider 23 receives the 
output signals from the low pass 13 and the low pass 22 and outputs a ratio of the two 
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input signals. The output of the divider 23 is connected to another low pass 24, which, 
in turn is connected to the phase/gain compensator 17, 

The operating principle of the receiver of Fig. 1 and a method for 
operating the receiver of Fig. 1 according to an exemplary embodiment of the present 
5 invention will be described in the following. 

Phase imbalance estimation 

The method for phase imbalance estimation according to an exemplary 
embodiment is based on the observation that after the demodulator, the I and Q 
components affected by imbalance can be written as function of the i and q components 
1 0 without imbalance as shown in equation (4): 
I = kd 

Q = k(-isincp + qcoscp) (4) 

Equations (4) hold true if the phase of the I-branch of the demodulator is 
aligned with the phase of the i component of the received signal. 
1 5 (i,q) is assumed to be the transmitted signal for UMTS, i.e. the signal 

without including the effects due to the propagation through the multi-path fading 
channel. From the complex scrambling code properties and from the modulation 
scheme adopted the following equations (5) can be inferred. 
<i-q>= <q-i> = 0 

20 <ii> = <qq> = P/2 (5) 

Thus the correlation matrix of (i,q) is diagonal as already described 

above and the terms on the principal diagonal depend from the transmission power. 
From (4) and (5) it can be written: 

<LQ> = -k-<i 2 >sin<p 

25 <F> - k-<i 2 > (6) 

<I-Q> is the output signal of the low pass 12 and <P> is the output signal 
of the low pass 13. Both signals are input to the means for generating a ratio 14. 

As shown above, it can be inferred that: 
sin<p = -<IQ>/ <F> (7) 
30 where <I*Q> is defined as cross correlation of I and Q, <P> as quadrature 

of I, and <Q 2 > as quadrature or square of Q. 
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When propagating through the multi-fading environment, the ideal 
QPSK (quadrature phase shift keying) constellation at the transmitter is 
rotated/attenuated/delayed for each path. This has an influence on the cross correlation 
properties of (i,q) before the receiver. Furthermore, the signal becomes intrinsically not 
5 stationary and the statistics of the (i,q) samples get more complicated by the fact that 
there can be correlations between delayed versions of the same signal. It has been 
shown above that these effects have no influence and the residual presence in the cross 
correlation in the received components (i r , q r ) is only due to the phase imbalance. 

For compensation of the phase imbalance, the resulting signal is put into 
10 phase/gain compensator 17, which applies the following equations: 
ic = I 

q c = (Q + I-sin<p) / coscp (8) 
where: 

I 

cos^ = ((</ 2 > 2 -<I>Q> 2 )l<I 2 > 2 ) 2 (9) 
15 ic and q c are the output signals of the phase/gain compensator 17 on the 

lines 18 and 19. 

Gain imbalance estimation in the presence of phase imbalance 
In the following, the method is extended to the case in which there is also 
an amplitude mismatch besides the phase imbalance. 
20 In case there is also an amplitude mismatch, the above equations can be 

written as follows: 
I = g(t)(M) 

Q = g(t)(k 2 (-isuKp + qcos<p) 

<I,Q> = ki • k 2 • <g(t) 2 > <q 2 >sin(cp) 
25 < IQ >= -fc, -k 2 < g(t) 2 >< i 2 > sin(?)) = -k x • k 2 < g(t) 2 >< q 2 > sin(p) 

<Q. Q> = k 2 2 • <g(t) 2 > (<i 2 >cos(<p) 2 + <q 2 >sin(q>) 2 )= k 2 • <q 2 > 

<[.I> = k 1 2 -<g(t) 2 ><i 2 > 

From these equations it follows: 

<Q> / (<sh^(£>) xl2 = -sin(<p) 
30 <Q 2 >/<I 2 > = k 2 2 /k l 2 (10) 
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Where the ratio <Q*> / <I 2 > is the output signal of the divider 23. After a 
low pass filtering by means of low pass 23, the ratio<Q 2 > / <^> is provided to the 
phase/gain compensator 17. 

Note that the temporal coefficient g(t) representing the amplitude of the 
5 channel dependence on time can be simplified, because it appears both in the numerator 
and in the denominator of equation (10), which can be considered a generalisation of 
equation (7). The phase factors are removed also, because they would anyway disappear 
in the auto-correlation matrix. 

Feed-back scheme for phase imbalance 
10 Fig. 2 shows a feed-back phase imbalance estimation and compensation 

device according to an aspect of the present invention. As shown in Fig. 2, an incoming 
complex signal r(t) is received at a receiver station not shown in Fig. 2 and is further 
processed such that the complex signal r(t) is provided to input connections 1 and 2 of 
the demodulator 3. After that, the complex signal r(t) is demodulated. Between the input 
15 connections 1 and 2, A cos(o)t) represents the waveform used in the receiver to 
demodulate the incoming signal. 

The I component on line 4 and the Q component on line 5 of the 
demodulated complex signal r(t), which are affected by phase imbalance or gain 
imbalance, are put through phase/gain compensator 17 and afterwards provided to lines 
20 18 and 19. The output signals of the phase/gain compensator 17 are designated with i c 
and q c . The ic and q c signals are fed into means 8 for taking the cross correlation of the I 
and Q components, such as a multiplier by means of lines 6 and 7. Additionally, the ic 
signal is fed into means 9 for taking the quadrature or square of the I component, such 
as a multiplier which multiplies the ic with itself. After that, the resulting signals are 
25 provided to low pass filters 12, 13 via lines 10, 1 1, respectively. 

For estimation of the phase imbalance, a ratio between the cross 
correlation of I and Q and the quadrature of I is generated by means for generating a 
ratio 14, such as a divider. After that, the resulting signal is integrated by means of an 
integrator 25 provided to line 16 and fed back into the phase/gain compensator 17 
30 located between lines 4, 5 and 18, 19, for compensating/estimating the phase imbalance 
or the gain imbalance of the I and Q components. Thereupon, the estimated and 
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compensated I and Q components i c and q c , respectively, are fed back to lines 18, 19 to 
be fed into means for synchronisation 20 such as a UMTS synchronizer. 

In the following, an exemplary embodiment of a method of the present 
invention for operating the device of Fig, 2 and its operating principle will be described. 
5 In the feed-back phase imbalance estimation and compensation device 

according to the present invention and the method for operating the feed-back phase 
imbalance estimation and compensation device according to the present invention, the 
compensated output of the phase/gain compensator 17 feeds the imbalance estimation. 
Some simplifications have been introduced in the formulas for the compensation and for 
10 the phase estimation. This is possible because in the next iteration in the loop the phase 
error gets smaller and the approximation of the sinus with the value of the angle and the 
cosines with 1 is appropriate. The principle of the method will be described and 
explained in the following: 

The following equation can be used for the estimation of the phase <p: 
15 q> = -<[-Q>/<I 2 > (11) 

The phase compensation carried out in the phase/gain compensator 17 
can be described with the following equation: 

ic = I (12) 

q c = (Q + I-<p) 03) 
20 A further simplification which can be introduced for the phase 

compensation feedback method is to use the following expression for the phase: 
<p = -k*<IQ> (un-normalized expression for the error) 

where k is a constant. This is possible because this term is equal to zero when the phase 
imbalance is zero. Hence, this the signal represented by this term can be then used as 
25 error in the feedback scheme. 

According to an aspect of the present invention, a more efficient way to 
implement the feedback scheme is via an iterative procedure. According to this aspect, 
the low pass filters and the integral may unified in a single computation step for 
implementation. This will be described with more detail in the following: 
30 Iterative implementation for phase imbalance 

In the following, i e (n), q e (ri) are the phase compensated received in- 
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phase components and quadrature components. 

1) At an instant n, the following error function is computed: 
e =ic(»Vtfc(*) 

2) Then, the following equation is used to ohtain a filtered integrated 
5 expression for the error: 

e f =e + p-e 

where the parameter p is chosen by matching convergence speed and stability. 

3) Then, the next (subsequent) i c (n + 1), q c (n + 1) is computed from the next 
i(n + 1), g(n + 1) received samples as follows: 

/ („ + l) = z(n + l) 

10 

g c (n + 1) = g(n + 1) + i(n + 1) • e f 

Then, the method returns iteratively to step 1). It has to be noted that the 

final value for ef represent the estimated sinus of the phase imbalance multiplied for the 

gain factor k. 

Feedback scheme for gain imbalance 
15 Fig. 3 shows a device for gain compensation with a feedback scheme. 

The signal ic on the line 18 is provided to means 30 for taking the 
quadrature or square via line 7. The output signal of the means 30 for taking the 
quadrature is filtered by means of a low pass 31 and then provided to means for 
subtracting 32 such as a subtracter. This signal is designated as signal A in Fig. 3. 
20 The signal q c on the line 19 is provided to means 9 for taking the 

quadrature or square via line 6. The output signal of the means 13 for taking the 
quadrature is filtered by means of a low pass 3 1 and then provided to the means for 
subtracting 32. This signal is designated as signal B in Fig. 3. The means for subtracting 
32 are configured such that the output signal of the low pass 31 is subtracted from the 
25 output signal of the low pass which is indicated by B - A in Fig. 3. 

The device depicted in Fig. 3 operates as follows: The received input is 
demodulated, and then gain compensated in the each iteration of the loop. The gain 
compensated components are squared, filtered and subtracted to build an error function 
which is integrated and then provided to the gain compensation block 17, to be 
3 0 processed together with the incoming signal. 

The equations that give a rationale of the method are the followings. 
Starting from the equations: 
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<IQ> 



(14) 
(15) 



(17) 



V</ 2 ><Q 2 > 

<Q 2 > = k 2 2 
<I 2 > k x 2 

wherein x can be defined as follows: 

k 2 -k\ <! 2 >-<Q 2 > ^ l-K'jk.-kJ 
* k 2 <I 2 > k 2 k, 

which becomes: 

. _ k 2 -k 2 <I 2 >-<Q 2 > _ 2-k r (Jc,-k 2 ) 2-{k x -k 2 ) 
k 2 <I 2 > " k 2 " k r 

and 

gc =\-Q.5.x = ^- (16) 
k 2 

where the following approximation has been used: 
10 The gain compensated components are obtained as follows: 

Note that the derivation is based on an approximation becoming more 
precise with a decreasing difference between k x and k 2 under the condition that it is 

realised in the next iterations of the loop. 
15 A further simplification is obtained by considering: 

x l = k x 2 - k\ (un-normalised expression for the error) 

for applying the method. In fact, this last expression is zero for gain balanced 
components, and does not involve the ratio operation which can present some problems 
in the presence of noise. 
20 Iterative implementation for gain imbalance estimation and 

compensation 

An even more efficient way to perform the low pass filtering and the 
integration operation is represented by the following steps in accordance with an 
exemplary embodiment of the present invention: 
25 1) In an instant n, the following term is computed: 

e =i e (*) 2 -<lc(n) 2 
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2) Then, the following equation is used to obtain a filtered integrated 

expression for the error: 

e f = e + p • e 

where the parameter p is chosen by matching convergence speed and stability 
5 3) Then, the gain g cab be derived as follows (k is a given constant): 

g = l-k-e f 

4) The computed gain is then applied to the next i and q samples: 

/ c (n + l) = /(n + l) 
? c (n + l) = ?(/* + l) g 

5) Then, the method iteratively restarts from step 1) 

1 0 Note the final value of g will represent the estimated value for the gain. 

Advantageously, this iterative method is valid in presence of a phase imbalance and in 
absence of phase imbalance. 

Phase and gain imbalance: feedback/iterative method. 

The two loops represented in Figs. 2 and 3 can be serially concatenated 

15 to obtain simultaneous compensation of phase and gain. In particular, the output of the 
loop represented in figure 2 can be the input of the loop represented in figure 3 and the 
final output would be compensated both in phase and in gain. This is a consequence of 
the fact that in the equations (10) phase and gain imbalances are independent. Also, if 
this independence does not apply in the case in which un-normalised expressions are 

20 used for the errors to reduce complexity, it can still be shown that the serial 

concatenation of the loops operates very efficiently. In fact, in this last case, the residual 
phase imbalance influence the estimation of the gain imbalance in the next iteration of 
the loops, but this influence becomes zero when the phase error becomes zero. 
Accordingly, the coupling created by using the un-normalised forms for the errors 

25 creates a reciprocal influence in terms of convergence speed, but does not affect the 
convergence itself. 

Figs. 1 1 and 12 depict results achieved with two concatenated iterative 
loops in which the phase and the gain are simultaneously compensated. 
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As will be shown in the following, the above iterative approach based on 
un-normalised errors allows for a very efficient and more general estimation and/or 
compensation of gain/phase imbalances. In fact, starting from the following equations: 

I = k x i (18) 

Q = -k 3 • i • sin <p 4* k 2 q • cos#> 

in which fc 3 is used for the part of i injected on the q branch, it can be 

shown that: 

< / . Q >= k x k 3 sin(p) < i 2 > 

<I 2 >-<Q 2 >=< i 2 > (k 2 - k 2 sen 2 <p) - fc 2 ? 2 cos? (19) 

As can be seen from these equations (19), the first one holds true if one 
of the two phase and gain unbalanced components is multiplied with an arbitrary factor. 
Also, <IQ> gives an error that is zero only if the phase imbalance is zero. On the other 
hand, in absence of phase imbalance, the second equation is zero, if and only if: 

<I 2 >-<Q 2 >=(<z 2 >k 2 -<g 2 >* 2 ) = 0; (20) 
<i 2 > _k 2 



<q > k : 



2 



In this way, a cascade scheme can be designed, in which the equation for 
i the phase is applied on phase unbalanced components as given in equation (1 1). The 
output of this first loop feeds a second loop, which applies the procedure given in 
equations (12) and (13) to estimate and compensate for the gain imbalance. This scheme 
will compensate both phase and gain imbalance in a computational effective way 
whenever they can be expressed in the form given by the equations (18). 
) Fig. 4 shows components of a UMTS receiver. There is a radio frequency 

(rf) demodulator 40 which output signal is input to an A/D converter 41. The output of 
A/D converter 41 is input to an RRC filter 42 which is connected to a down sampler 43. 
From the down sampler 43, the signal is output to a further processing. The lines 44, 45 
and 46 indicate positions within the UMTS receiver, where the above methods for phase 
5 and/or gain imbalance estimation and/or compensation may be applied. As can be seen 
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from Fig. 4, phase and/or gain imbalance estimation and/or compensation may be 
applied between the A/D converter 41 and the RRC filter 42, i.e. before the RRC filter 

42, between the RRC filter 42 and the down sampler 43, i.e. before the down sampler 

43, and after the down sampler 43 before the further processing. Therefore, according to 
5 the present invention, 1 a great degree of freedom is provided for the design of a receiver 

according to the present invention. 

Performance of the inventive method and device 

Figs. 5-7 show a ratio between the estimation error and the true value for 
filters of different length as resulting from a simulation based on typical conditions for 

10 UMTS receivers extracted from the standard specifications. The phase imbalance error 
reaching from 0 to 5 degrees and the length of the low pass 2 in Fig. 1 are represented in 
the X, Y axes. The unit for the length of the low pass filter is 1 pilot symbol, i.e. 256 
chips. The length of the low pass filter 1 in Fig. 1 is always assumed to be 256 chips. 
The Z-value E represents the square root of the variance of the estimated value of the 

15 actual imbalance. It is related to the mean magnitude of the residual phase error after 
compensation. 

Li Fig. 5, the conditions considered correspond to 4 paths, a velocity of 3 
km/h and lorx/loc = -3 dB. 

In Fig. 6, the conditions considered correspond to 4 paths, a velocity of 
20 120 km/h and lorx/loc = -3 dB. 

In Fig. 7, the conditions considered correspond to 4 paths, a velocity of 
250 km/h and lorx/loc = -3 dB. 

All three simulations show very good results for E for a filter length of 
low pass filter 2 in Fig. 1 of more than 50 pilot symbols or a phase error of more than 2 
25 degrees. 

Fig. 8 shows the effect of the phase imbalance on the BER (Bit-Error 
Rate) in a case in which AWGN propagation conditions (Additivd White Gaussian 
Noise) are considered and no other losses are included. Between 0 and 4 degrees the 
BER increases rather slowly. From 4 to 7 degrees, the increase of the BER is moderate 
30 and from 7 to 10 degrees, the BER increases fast. As can be seen in the range 0-5 
degrees, the phase imbalance can imply an increase in the BER of about 60%. 
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Fig. 9 shows the performance which can be achieved in typical 
conditions. The phase error decreases exponentially with the integration length. By 
choosing an integration length of 16 x 256 chips, corresponding to 16 pilot symbols, the 
phase error is reduced in mean of the 90%. 
5 Fig. 10 shows results related to a case in which the phase is changing at a 

fixed frequency of 5 Hz. This case is used to give a representation of the tracking 
capability of the method. The conditions considered correspond to 4 paths and a 
velocity of 120 km/h. The effect of the noise is included, inter-cells interference at + 3 
dB with respect to the desired signal. The integration period is 16 slots and the mean 

10 square error achieved is 0.57 on a peak phase error of 5 degrees. 

Figs. 11 and 12 show the performances of the iterative procedure applied 
when both phase and gain imbalance are present. Fig. 1 1 shows an estimated phase with 
speed =3 km/h and 3 equivalent paths in the iterative application of the method at the 
location indicated with arrow 44 in Fig. 4. The initial condition for the gain imbalance 

15 is 0.8 and 5 degrees for the phase. The term k 3 in the equations (18) is equal to 0.5. 

Figure 1 shows the error for the gain with speed =3 km/h and 3 
equivalent paths in the iterative application of the method at the location indicated with 
arrow 44 in Fig. 4.The initial condition for the gain imbalance is 0.8 and 5 degrees for 
the phase. The term fc 3 in the equations (18) is equal to 0.5 

20 As can be taken from Figs. 1 1 and 12, after approximately 300 pilot 

symbols, i.e. 2 frames corresponding to an interval of 20 msec, errors of 5 degrees in 
phase and 0.8 in gain are reduced to about 0.02 degrees for the phase and 0.01 for the 
gain. 

The method and device of the present invention is particularly well suited for 
25 UMTS/W CDMA, where a complex scrambling code is used to modulate the information 
symbols. However, the method and device of the present invention may also be applied for 
any transmission with QPSK modulation, a modulation scheme based on a complex 
scrambling code and a frequency selectivity of the channel with respect to transmission 
bandwidth. 

30 The device according to the present invention may be designed by means of 

EPLDs , such as the Altera EPLDs or by means of suitable calculation devices. The 



WO 2004/054194 ^^T/IB 2003/005 779 

20 

methods of the present invention may be implemented as software program code in a 
suitable programming language such as C++ and stored on a suitable computer readable 
storage device such as a compact disc. 



